Fourth-order perturbative extension of the singles-doubles coupled-cluster method 
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Fourth-order many-body corrections to matrix elements for atoms with one valence electron are 
derived. The obtained diagrams are classified using coupled-cluster-inspired separation into con- 
tributions from n-particle excitations from the lowest-order wavefunction. The complete set of 
fourth-order diagrams involves only connected single, double, and triple excitations and discon- 
nected quadruple excitations. Approximately half of the fourth-order diagrams are not accounted 
for by the popular coupled-cluster method truncated at single and double excitations (CCSD). Ex- 
plicit formulae are tabulated for the entire set of fourth-order diagrams missed by the CCSD method 
and its linearized version, i.e. contributions from connected triple and disconnected quadruple ex- 
citations. A partial summation scheme of the derived fourth-order contributions to all orders of 
perturbation theory is proposed. 

PACS numbers: 31.15.Md,31.15.Dv,31.25.-v 
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I. INTRODUCTION 

Atomic tests of the low-energy electroweak sector of 
the standard model require both high-precision measure- 
ments and ab initio calculations of matching accuracy. 
The most precise measurement to date of parity violation 
in atoms has been carried by Wieman and co-workers us- 
ing i^'^Cs. The accuracy of this experiment is about 
0.4%, while the relevant theoretical quantity is calculated 
with 0.4-1% uncertainty, depending on the authors's esti- 
mates ^ . A keen interest in reducing the uncertainties 
is stimulated by a possible deviation of the resulting nu- 
clear weak charge from the prediction of the standard 
model. This deviation was first reported in Ref. and 
then scrutinized in Refs. [||. Very recent analyses [|| of 
parity violation in ^'^^Cs focused on effects of the Breit in- 
teraction, vacuum polarization, and neutron "skin" , each 
contributing at the level of 0.2-0.6%. However, the ef- 
fects of higher-order correlations beyond those consid- 
ered in high-precision calculations by Dzuba et al. ^ 
and Blundell et al. j|] remain to be understood. Here 
we discuss in detail a possible extension to the method 
employed in Ref. j^. 

The key to the 1% accuracy achieved in Refs. ||] 
lies in the application of all-order methods based on rela- 
tivistic many-body perturbation theory (MBPT). These 
techniques, although summing certain classes of MBPT 
diagrams to all orders of perturbation theory, still do 
not account for an infinite number of residual diagrams. 
It seems natural to augment a given all-order technique 
with some of the omitted diagrams so that the formalism 
is complete through a certain order of MBPT. To illus- 
trate, the random-phase approximation (RPA) fully 
recovers second order matrix elements but does not sub- 
sume all third-order diagrams. Among the omitted third- 
order contributions so called Brueckner-orbital diagrams 



are known to be numerically as important as the RPA 
sequence (see, e.g., discussions in Refs. 0, ||)- 

By the same virtue, certain diagrams starting from 
the fourth order of MBPT are missed in the popular 
coupled-cluster expansion ||, |o[ |ll| truncated at the sin- 
gle and double level of excitations (CCSD), although all 
third-order contributions are recovered [|l2|. It has been 
shown that one of the subsets of the fourth-order 
terms missed by the CCSD method does contribute as 
much as a few per cent to Cs hyperfine-structure con- 
stants. At the same time, the considered subset leads 
to worse theory-experiment agreement for electric-dipole 
amplitudes [Q. We anticipate that a systematic ac- 
counting of all omitted fourth-order contributions to ma- 
trix elements in the CCSD method may lead to more 
accurate ab initio results. Here we derive such compli- 
mentary fourth-order many-body contributions for ma- 
trix elements. 



The paper is organized as follows. Basic starting 
formulas and notation of many-body perturbation the- 
ory (MBPT) are introduced in Section II. The linked- 
diagram expansion specialized to atoms with a single va- 
lence electron is discussed in Section [II. The derived 
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wave-functions through the third order of MBPT are 
discussed in Section IV and their relation to the trun- 
cated coupled-cluster method in Section Finally, the 
derived fourth-order corrections to matrix elements are 
tabulated in the Appendix and classification of the di- 
agrams is given in Section VI. Fig. || summarizes the 
results of our work. 



The fourth-order expressions presented here may be 
useful for an analysis of completeness of all-order meth- 
ods and for designs of a hierarchy of next-generation ap- 
proximations in atomic many-body calculations. As an 
example, we discuss all-order generalizations of the de- 
rived fourth-order contributions. 
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II. PARTITIONING OF THE ATOMIC 
HAMILTONIAN 

Here we briefly recap starting formulas of many-body 
perturbation theory (MBPT) for atoms with one valence 
electron. Our derivation in the fourth order of many- 
body perturbation theory may be considered as an exten- 
sion of the work by Blundell et al. They presented 
formulas from first-, second-, and third-order perturba- 
tion theory. For the convenience of the reader we keep 
most of the original notation from Ref . . 

The many-body Hamiltonian of an atomic system may 
be represented as 



^ /inuc(ri) -I- ^ J7HF(ri 

\ i i 



(1) 



where /inuc includes the kinetic energy of an electron 
and its interaction with the nucleus, Uhf is the Hartree- 
Fock potential, and the last term represents the residual 
Coulomb interaction between electrons. The summations 
go over all electrons in the system. In MBPT the first 
part of the Hamiltonian is treated as the lowest-order 
Hamiltonian Hq and the residual Coulomb interaction as 
a perturbation V7 . 

For atoms with one valence electron outside a closed- 
shell core the many-body wavefunction in the lowest or- 
der is a Slater determinant constructed from core 
and valence single-particle orbitals Uk which satisfy 



(/inuc(r) + Unpir)) Ufc(r) = efcUfe(r) 



(2) 



The solutions of the above one-particle equation form a 
basis for application of the formalism of second quantiza- 
tion. In the second quantization the lowest-order Hamil- 
tonian Ho and the perturbing residual Coulomb interac- 
tion Vj may be expressed as 

Hq = '^Eialai , (3) 

i 

ijkl ij 

where a| and Oi are creation and annihilation operators 
for a one-particle state i. 

The Coulomb integral gijki is conventionally defined as 

g^Jkl = I ul (r) u] (r') y^^uk (r) m (r') d\ d\' . (5) 

The matrix elements of the Hartree-Fock potential may 
be expressed in terms of the antisymmetrized Coulomb 
integral (jijki = giju - gijik as 



9iaja 



(6) 



Here the summation is over core orbitals; this potential 
is the so-called frozen-core Hartree-Fock potential, i.e., 
first the core orbitals are calculated employing the self- 
consistent Hartree-Fock procedure and then the rest of 
the one-particle states are obtained using Eq. (^) without 
varying the determined core orbitals. Finally, in the lan- 
guage of second quantization the lowest-order wavefunc- 
tion corresponds to = a!j^|Oc), where v labels the 
one-particle state of the valence electron and the quasi- 
vacuum state |0c) describes the closed-shell core. 

^From a practical standpoint derivation of MBPT ex- 
pressions is greatly simplified by the introduction of nor- 
mal form of the operator products N[. . .] and by a subse- 
quent application of the Wick theorem |ll|. The notion 
of normal products arises from separation of one-particle 
states into two general categories - occupied in the quasi- 
vacuum state |0c) ( i.e., core orbitals enumerated by let- 
ters a,b,c,d) and complementary excited states (indices 
m,n,r,s). Unspecified orbitals are labelled by indices 
i,j,k, and I. In this scheme the one-particle valence 
states V and w are classified as excited orbitals. 

With the normal products 



and 



(7) 



(8) 



ijkl 



where e'^c'' = Eo^^ and e'^c''^ = -\Y.a (^HF)oa; in the 
following discussion we omit these nonessential offset con- 
tributions. 

It is worth noting that there is no one-body part of 
the perturbation Vi present in Eq. (H); this fact demon- 
strates the utility of the frozen-core Hartree-Fock poten- 
tial in MBPT. In Ref. the case of a model poten- 
tial differing from J7hf was investigated explicitly and it 
was found that the number of resulting diagrams is sub- 
stantially larger than in the Hartree-Fock case. Due to 
the very large number of diagrams in the fourth order, 
here we restrict our attention to the practically impor- 
tant frozen-core Hartree-Fock case. 



III. LINKED-DIAGRAM EXPANSION 

We proceed to the derivation of many-body contri- 
butions to wavefunctions using the formalism of the 
generalized Bloch equation The Bloch equation 

is formulated for the wave operator which relates 
the exact wavefunction \^v) to the lowest-order result 

ai\Qc) as 



1*1°)) 



(9) 



It should be noted that as defined, this exact wavefunc- 
tion is not normalized, rather an intermediate normaliza- 
tion scheme (^i"'!*!,) = 1 is employed in the formalism. 
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The exact correlation energy of the one-valence electron 
system is given by 



(0)\ 



(10) 



The wave-operator satisfies the linked-diagram version 
of the generalized Bloch equation 



i^v, Ho] = {Q Vi - {n, -i)PVi rijii 



linked ' 



(11) 



where the operator P — 1^1, | projects on the 

lowest-order wavefunction and Q = 1 — P is a com- 
plementary projection operator. The subscript "linked" 
in the above equation prescribes that all the unlinked 
Brueckner-Goldstone diagrams are to be discarded; a di- 
agram is said to be unlinked if it contains a disconnected 
part with no free lines other than valence lines. Finally, 
[f2„, Hq] is the commutator Qy Hq — Hq fly. 

The traditional Rayleigh-Schrodinger perturbation 
theory is recovered from the Bloch equation (^l]) by ex- 
panding the wave operator in powers of the residual in- 
teraction Vi, fly = J2n=o ^t'"''- The resulting recursive 
relation is ITl|l 



m=l J 



,(12) 



linked 



Here the iterations start with fl^^ = 1. A corresponding 
perturbative expansion of correlation energy reads 



6Ey = Y.5E^-^ = ^(vI/(«)|Fzl7("-i)|vI/(' 



(0)\ 



(13) 



n=l 



The last term on the r.h.s. of Eq.(|iJ) gives rise to 
so-called "folded" or "backward" diagrams ||ll|. Instead 
of calculating the explicit contributions of folded dia- 
grams we use an all-order approach which incorporates 
their effect in modified energy denominators. Such a 
reformulation allows for a direct link to the coupled- 
cluster method outlined in Section |^. The exact wave- 
operator r2„ may be separated into valence and core 
parts, ^y = ni""^ + VL"""", the ^l""^ part promoting a va- 
lence electron from the j^i"'') determinant into an excited 
state, ri'^"™, describing excitations of core electrons, does 
not depend on any particular valence state. Similarly, the 
correlation contribution to the total energy of the system 
5Ey may be broken into corrections to the energies of va- 
lence and core electrons, 5Ey — SEy^^ + SE'^™'^. 

Suppose that the valence removal energy £„ -I- SEy'^^ is 
known at the desired order of perturbation theory (e.g., 
from coupled-cluster calculations) or from experiment. 

Projecting the original Bloch equation ( pT| ) onto l^'i"') 
and using the definition of the projection operator P to- 
gether with Eq. (|l^) for the correlation energy, one may 
show that 



{QViny},,^ 



kcd 



Notice that the last term is represented by a product of 
two valence contributions, since all other terms produce 
unlinked diagrams. Expanding the commutator and ex- 
plicitly breaking the term {Q ^/^^ulunkcd '^^^'^ valence and 
core contributions we arrive at 

{sy + 5E:-' - Ho) f^r'l*i°^> - {{Q Vi f^Jnnkcd)''' 1*^'"^): 

{Sy - Ho) 17™-!*!°)) = {{Q Vl f7.}unked)'°" ■ 

Accounting for the folded diagrams in this way leads 
to an additional shift JE'™' in energy denominators of 
diagrams for the valence part of the wave operator fi™'. 
Mnemonically, every occurrence of the Hartree-Fock en- 
ergy of the valence electron Ey in the energy denominators 
has to be replaced by the total removal energy -l-Ji?™', 
since (e„ - Ho) 17™'°!*!°)) simplifies to (-Hq) n''°"'\Oc). 
Keeping this rule in mind, we may combine the above 
equations 



{Sy + {5E:'^')-Ho) ny\¥^^) = 



{QViny} 



linked 



1*1°)) 



(14) 



where means that the (5i?™' correction should be 

included for the valence diagrams of i^y and discarded 
otherwise. 

We expand the wave operator in powers of the residual 
electron-electron interaction V7, Qy — X]n=o ^^"^ ^^'^ 
obtain 

(sy + iSE^') - ffo)f7l"+'^|*i"^) = 

^ linked 

with ril*^) = 1. This equation may be interpreted as a 
linked-diagram version of the Brillouin-Wigner perturba- 
tion theory for atoms with one valence electron outside 
a closed core. Introducing the resolvent operator 



Ry = {Ho - [ey + {6E^^')]) ' , 
we obtain (with |*1")) = 1]1")|*1°))) 

{g^/l*l"-^^)}^ 



(15) 



* 



(")\ _ 



linked 



(-l)"({i?.Qyz}unked)"l*i"^> 



^From this recursion relation we may generate correc- 
tions to wave functions at any given order of perturbation 
theory. In practice, the derivation is rather tedious and 
error-prone. We employed the symbolic algebra system 
Mathematica to derive the expressions presented in 
this work. 



IV. WAVEFUNCTIONS THROUGH THE 
THIRD ORDER OF MBPT 

For the derivation of fourth-order matrix elements one 
requires contributions to wavefunctions through the third 
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order. Expressions for = f^t,"'' through the 

second order may be found in Ref . [Esl . Although we fuUy 

derived !'!'«), to keep the manuscript to a manageable 
size, we present below only a qualitative discussion of the 
third-order correction to the wave-function. 

The contributions to the wave operator fly are con- 
ventionally classified by the number of excitations from 
a reference determinant = aj|0c). The first-order 

result, contains only double excitations drawn in 

Fig. § We may distinguish between valence and core 
excitations. The former promote the valence electron to 
an excited state (5^™') and the latter do not modify the 
state of valence electron (17'^°''°). With such a classifi- 
cation the diagram Fig. |l|(a) represents core doubles Dc 
and Fig. ^(b) valence doubles Dy. 



(a) Dc 



(b)D, 



FIG. 1: Brueckner-Goldstone diagrams for the first-order 
wave operator Horizontal dashed lines represent resid- 

ual Coulomb interaction between electrons and vertical lines 
are particle/hole lines. The valence line is marked by double 



o:V V::V Ql V: 



(a) S, 



(b) D, 



(c) 



(d) D, 



v.: 



V 



V.V 



(e) Tc 



(f) T. 



(g) Dnl 



FIG. 2; Sample contributions to the second-order wave op- 

(2) 

erator f2„ . 



We classify the disconnected quadruple contribution |^(g) 
as a nonlinear contribution of double excitations to wave- 
functions. 



V.V 




N[ 



® 



FIG. 3: A sum of the the quadruple contribution Fig. g(g) 
and a similar diagram with the order of two interactions re- 
versed factorizes into a normal product of double excitations. 
On the r.h.s. the energy denominators are to be evaluated 
separately. 



Several contributions to the third-order wave-operator 
fl^^ are shown in Fig. IJ. Single and double excitations 
shown in Fig. H(a-e) contain intermediate triple excita- 
tions. Diagraniy(f) is due to intermediate second-order 
quadruple excitation. For the sake of comparison with 
the coupled-cluster method we classify diagram ^(a) as 
the effect of core triples on core singles (S'c[Tc]), (b) as 
modification of core doubles by core triples {Dc[Tc]), (c) 
as the effect of core triples on valence doubles {Dy[Tc]), 
and (d) and (e) as the effect of valence triples on valence 
singles and doubles (S'i,[71,] and Dy[Ty]). Finally, dia- 
gram ^(f ) may be classified as an effect of nonlinear dou- 
bles, Fig. ||(g), on valence doubles {Dy[Dni])- It is worth 
noting that the third-order wavefunction contains con- 
nected quadruple excitations and some additional discon- 
nected excitations; these corrections do not contribute to 
the fourth-order matrix elements. 





(a) ScCTc] (b) Dcffc] (c) Dvffc] (d)Sv[rv] 




0' 



V. 



(e) DvCTv] 



(f) DvPnl 



(g) 



(2) 

The second-order operator ilj, contains excitations up 

(21 

to quadruples. Examples of contributions to Q V are 
drawn in Fig. ||. Diagrams ^(a) and (b) represent some 
of the second-order core singles and doubles. Valence sin- 
gles and doubles are drawn in Fig. ||(c) and (d) respec- 
tively. Diagrams ||(e) and (f ) represent core and valence 
triple excitations, and (g) — disconnected quadruple ex- 
citations. A sum of the the quadruple contribution ||(g) 
and a similar diagram with the order of the two interac- 
tions reversed is known to factorize into a normal product 
of double excitations [jll] ; this is demonstrated in Fig. ||. 



FIG. 4: Representative diagrams for the third-order wave 
operator Sl^'''. 



V. COUPLED-CLUSTER METHOD 

The coupled-cluster (CC) formalism H, |l^ is widely 
employed in atomic and nuclear physics, and quantum 
chemistry [ |T7[ . The main goal of the present work is 
to identify fourth-order contributions to matrix elements 
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not included in the truncated singles-doubles coupled- 
cluster method, and here we review the relevant features 
of this all-order approach. 

The key point of the coupled-cluster method is the in- 
troduction of an exponential ansatz for the wave opera- 
tor ITTi 



n = iV[exp(5)] = 1 + S+ ^N[S^] 



(16) 



The cluster operator S — ilconn is expressed in terms of 
connected diagrams of the wave operator ft, an example 
of disconnected diagram being Fig. ^(g) . The operator S 
is naturally broken into cluster operators Sn combining n 
simultaneous excitations from the reference state j^'i"') 
in all orders of perturbation theory. 

Let us specialize the general formalism of Ref. to 
the case of atoms with one valence electron. A set of 
coupled equations for the cluster operators may be found 
by considering connected diagrams on both sides of the 
modified Bloch equation dll) 



(e„ -I- (JK'O - ^o) Sn^iQ Vi n,}^ 



(17) 



where 6Ey is determined by Eq. iToh and wave opera- 
tor by Eq. Term (^£;™'jaccounts for folded 
diagrams; it is to be omitted for core and included for 
valence clusters. Successive iterations of such all-order 
equations explicitly recover order-by-order MBPT con- 
tributions to the wave operator discussed in the previous 
sections. 

In most applications the full operator S is truncated 
at single and double excitations (CCSD method). For 
univalent atoms the CCSD parameterization may be rep- 
resented as 



S' 



SD 



Si + S2 — 

Et ^ t t 

Pnia ^rn^o, ^ ^ / ^ Pmnab ^ra^n^b^a ~t- 

ma mnafc 

^ ^ Pmv ^rri^v ^ ^ Pmnva ^^^n^o,^v ; (-1^) 
m^v mna 



where the first two terms represent single and double 
excitations of core electrons and the remaining contribu- 
tions are valence singles and doubles. 

It is worth emphasizing that the CCSD method is 
an all-order method. For example, first-, second- 
and third-order diagrams Fig. |^(b). Fig. ^(d), and 
Fig. ^(g) are encapsulated in the valence doubles term 
Smna Pmnva al^al^aaay . Similarly, the CCSD method ac- 
counts for all single and double excitations (both core 
and valence) shown in Figs. ^ and ^. At the same 
time connected triple and higher-rank excitations are 
not accounted for by the CCSD method, examples be- 
ing Fig. |2|(e),(f) and Fig. ^ (a~e). Although diagrams 
Fig. ^ (a-e) are nominally single or double excitations, 
they contain connected triples as intermediate excitations 
and are not included in the sequence of CCSD diagrams. 



A Zmearized version of the CCSD method (LCCSD) is a 
further simplification of a hierarchy of all-order methods 
based on the coupled-cluster formalism. In this approx- 
imation SlJ;^'^^^ = 1-1- 5^°. For alkali- metal atoms the 
LCCSD method was employed in Refs. |l|, |l|, |l|, @. 
Compared to the full CCSD approximation the linearized 
version misses a subset of diagrams shown in Fig. |^(g) 
and 1(f). 

To reiterate, connected triple excitations and discon- 
nected quadruple excitations first appear in the second 
order wavefunctions. In order to systematically extend 
the CCSD method one has to investigate the contribu- 
tions of connected triple excitations and the role of non- 
linear contributions for the linearized CCSD approxima- 
tion. 



VI. MATRIX ELEMENTS 

We investigate the fourth-order corrections to matrix 
element of a one particle operator Z = z{ri). In sec- 
ond quantization 



(19) 



where N[. . .] denotes normal form of operator products. 
We are mainly interested in matrix elements of non-scalar 
operators, like electromagnetic transition amplitudes or 
pseudo-scalar operators, like the electroweak interaction. 
For such operators the contribution from the zero-body 
term Zaa vanishes and we disregard it in the following 
discussion. 

The exact matrix element between two valence states 
w and V is given by 



(20) 



where flu, and Qy correspond to wave operators for 
valence states w and v respectively. Since the wave- 
operators were derived using the intermediate normal- 
ization scheme, we introduced normalization factors 

iV. = (vI/(°)|f]tf^„|vI/W) 

in the definition of matrix element. 

BlundcU et al. have demonstrated that discon- 
nected diagrams in the perturbative expansion of the nu- 
merator and the denominator of Eq. ( 20|) cancel. Their 
final expression for the exact matrix element reads 



(^"'■")conn 



171 '(21) 



{[l+mcon„] [1 + W)conn]} 

where 

z™- EE {o,\nlzny\o,) = {o,\ (r!-'-)t zr!™-|o,) 
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and the remaining contributions of Z„,.„ = 
{Oc\aii]Ul^ Z ilyal\Oc) are grouped into the valence 
part The dia grams of exphcitly depend 

on valence indices w and v. The valence part of the 
normalization factor A^™' is defined in a similar fashion. 
The core contribution Z'^°'^'^ vanishes for non-scalar (and 
pseudo-scalar) operators and we disregard Z'^™° in the 
following discussion. Notice that all the diagrams in 
Eq. ( pT| ) must be rigorously connected as emphasized by 
subscripts "conn". 

The formulas for contributions to matrix elements 
through the third-order of MBPT were presented in 
Ref. The linearized coupled-cluster approach trun- 
cated at single and double excitations (LCCSD) fully re- 
covers the matrix elements through the third order [p^ . 
Here we investigate the contributions at the fourth order 
missed by the LCCSD method. 

To derive the fourth-order correction to a matrix ele- 
ment, we expand the matrix element and normalization 
factors into powers of the residual Coulomb interaction 



(fc) 



~2h^Q Nj}'\ Further, we employ 
the all-order result, Eq. (|2l|), and expand the normaliza- 
tion denominator into series. The result is 



WV, norm ' 



(*W|Z|vl/(3))- 



(*(? 1^1^(0) 



1^ val.c 



(22) 



where only connected valence contributions are to be 
kept. The normalization correction is given by 



wv, norm 



■^1,-1. v,r.vm rt \ V ' W } \ wv I 

Z \ / val,conn V / val.conn 

-^(^f^+^i^') , (23) 



val,conn 



where we used that Nij''' = and Z^v = z^v, 
element in the Hartree-Fock approximation. 

As we proceed to the derivation of the fourth-order 
diagrams we notice that the second line of Eq. (^) is 
the hermitian conjugate of the first line with a swap of 
valence indexes w and v. This observation allows us to 
consider only half of the diagrams since in numerical eval- 
uation the conjugated terms do not require additional 
programming efforts. 



the matrix 



VII. DISCUSSION OF FOURTH-ORDER 
DIAGRAMS 

We fully derived the fourth-order correction to ma- 
trix elements using Wick theorem. A set of simplifica- 
tion rules was implemented with the symbolic algebra 
system Mathematica [p^ . Excluding the normalization 
correction and folded diagrams, the resulting number of 
diagrams in the fourth order is 262. We counted both di- 
rect and all possible exchange forms of a given diagram 



as a single contribution. We excluded hermitian conju- 
gated terms from the counting procedure. The linearized 
coupled cluster approach, truncated at single and double 
excitations (LCCSD) recovers approximately half of the 
fourth-order contributions. The remaining diagrams are 
due to triple excitations (128 terms) and nonlinear contri- 
bution of double excitations (14 terms). Explicit expres- 
sions for these complementary contributions are given in 
the Appendix. 

We break all fourth-order contributions complemen- 
tary to the LCCSD subset of diagrams into nine classes: 



(Mifj) =^lx2(r«) + ^lx2(T, 

V /non-LCCSD 

Zox3{S4T,]) + ZoxsiD^T,]) + 
ZoxsiSciTc]) + ZoxsiDAT,]) + 

Zlx2{Dnl) + Zox3{Dnl) + ^norm(7t)) . 



(24) 



The representative diagrams for each class of contribu- 
tions are shown in Fig. ||. Here the diagrams Zix2i- ■ •) 
arise from evaluation of expression {¥^^\Z\¥y^) and its 
hermitian conjugate with a swap of valence labels w 
and V. Similarly Zqx3{- ■ ■) terms are generated from 
{¥°^\Z\-^i^^) + c.c. Finally (. . .) are due to nor- 

malization correction, Eq. (p3[). Further, we classify the 
diagrams by the presence of core triples (Tc) or valence 
triples {Ty). For Zqx3{- ■ ■) terms triple excitations oc- 
cur as an intermediate contribution (see Fig. |^) and we 
distinguish the effect of triples on lower-rank excitations, 
e.g. Dy[Tc] is the effect of core triples on valence dou- 
bles. Finally, the diagrams marked Dni are due to the 
effect of disconnected quadruple excitations. These di- 
agrams may be simplified to a direct product of double 
excitations, as demonstrated in Fig. ^j. 

The introduced classes of diagrams are illustrated in 
Fig. 1^. The numbers of contributions in each class are 
also given in that figure. Let us make some observations. 
First of all, none of the diagrams contain computationally 
intensive Coulomb integrals involving four particle states, 
e.g., Qmnrs- We also notice the absence of term Dc\T^, 
i.e., the effects of core triples on core double excitations. 
The core triples also do not contribute to the normaliza- 
tion correction. All these simplifications may lead to a 
design of an efficient numerical evaluation scheme. 

The dominant number of diagrams is due to valence 
triple excitations, the set Zix2{Tv) accounting for 44 
and the set ZQx3{,Dy\Ty]) for 36 contributions. We fur- 
ther distinguish second-order triples T by the nature of 
the orbital line connecting upper and lower interactions 
T = TP + T'',TP standing for a particle line and T'' for 
a hole line as illustrated in Fig. ^. Such a separation is 
motivated by considerations of computational complex- 
ity: the dia gram, involving summation over a small 
number of core states, may be calculated much faster 
than a similar contribution. We write 



Zlx2{Tv) 

ZoxsiDv [Tv]) 



Zix2iTP) + Zi, 
Zox3iD,[Tn)^ 



2(T^, 

Zox3iD4T;:]) . 
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Zi,2 (Tv) , 44 Zi,2 (Tc) , 20 Zox3 (Svffv] ) , 





Zox3 (Dv ffv] ) , 36 Zox3 (Dv [Tc] ) , 12 Z ox3 (Sc [Tc] ) , 



.0" 




® 



-o: 



Z 0x3 (D nl) J 7 Z 1x2 (D nl) ( Z norm (Tv ) / 8 

FIG. 5; Sample fourth-order diagrams involving triple ex- 
citations and non-linear coupled-cluster contributions. The 
one-particle matrix element is denoted by a wavy horizontal 
line. See the explanation in the text for diagram classifica- 
tion. The number of contributions for each class of diagrams 
is also shown; direct, all possible exchange, and the conju- 
gated graphs of a given diagram were counted as a single 
contribution. 



The formulas in the Appendix are grouped according to 
this scheme. 



v_. 



T ^ 



v_. 



FIG. 6: Separation of triple excitations based on a nature of 
an orbital line connecting upper and lower interactions. 
diagram involves particle line and a hole line. Similar 
separation may be carried out for core triples. 



as diagrams 
in Refs. [[| 



(3) 

The effect of triples on single excitations in VLy , such 
and (d), has been treated previously 

^ ^ Corresponding contributions to Mwv , 

Zoy^'i{Sv\Tv\) and Zo^^{Sc[Tc]) are shown in Fig. H It 
was found that this effect contributes as much as 5% to 
hyperfine-structure constants in Cs, and brings the ah 
initio calculations into 0.5% agreement with experiment. 
At the same time the experiment-theory agreement be- 
comes worse for electric-dipole matrix elements when the 
S\T] effect is included. To fully understand the role of 
triple excitations it is important to investigate all the 
enumerated effects on triple excitations, i.e. effect of 



triples on valence doubles, direct contribution of triple 
excitations to matrix elements entering Zix2, and also 
the normalization correction due to valence triple excita- 
tions. 

The linearized coupled-cluster method (LCCSD) 
[T^ , ^ [isll additionally disregards nonlinear terms in 
the coupled-cluster expansion. Therefore contributions 
to Sll^'' similar to one shown in Fig. ^(f) are omit- 
ted. These nonlinear terms lead to additional corrections 
Zox3(-D„;). A similar effect omitted in the LCCSD ap- 
proach is a direct contribution of disconnected double ex- 
citations to matrix elements represented by the diagrams 
of Zi-x_2{Dni) class. It is worth noting that considera- 
tion of the nonlinear contributions is key for accounting 
for the full set of random-phase-approximation diagrams 
with the CCSD method. 

We further notice that in the framework of tradi- 
tional Rayleigh-Schrodinger perturbation theory there 
are also contributions from so-called folded diagrams, as 
discussed in Section III. These folded diagrams origi- 



nate from the second-order valence energy correction (the 
first-order correction is zero in the frozen-core Hartree- 
Fock basis). Since both the CCSD method and its 
linearized version fully recover the second order ener- 
gies fl^ , in our approach we have omitted contributions 
of the folded diagrams. 

Finally, we would like to comment on a possible all- 
order extension of the derived fourth-order contributions. 
Ideally, the entire fourth order set of diagrams would 
be recovered by fully treating the triple and nonlinear 
double excitations within the traditional coupled-cluster 
approach. However, at the present state of computer 
technology such a full treatment hardly seems feasible in 
relativistic calculations. At the same time the coupled- 
cluster expansion truncated at the single and double ex- 
citations (CCSD method) presents an attractive start- 
ing point. The triple excitations may be treated semi- 
perturbatively, i.e., the triple excitations are replaced by 
a combination of "bare" Coulomb interaction and an all- 
order CCSD double excitation |p^ . 

The following modifications of the CCSD 
method should be made to partially sum the 
derived diagrams to all orders of perturbation 
theory: (i) Four classes of the derived diagrams 

may be accounted for by amending the traditional CCSD 
equations with a semi-perturbative contribution of triple 
excitations. Two of the desired modifications, S'^jTli] 
and S'ciT'c], were considered previously in Ref. jl3[ Q. 
(ii) In the diagrams 2'ix2(J\,), ^1x2(71;), and 2'„orm(T-o) 
the bottom and the upper (closing) Coulomb interac- 
tions should be replaced by all-order double excitation 
amplitudes. This generalization follows from consid- 
ering the relevant contributions in the coupled-cluster 
method, (iii) In the Zix2{Dni) diagrams all the Coulomb 
interactions should be replaced by all-order double ex- 
citation amplitudes, (iv) The linearized coupled-cluster 
expansions should include terms nonlinear in double ex- 
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citations to recover the diagrams ZoxaiDni) in all-order 
fashion. 



VIII. CONCLUSION 

An improvement of the accuracy of ab initio Coulomb- 
correlated calculations is necessitated by the latest ex- 
perimental and theoretical progress in studies of parity 
violation in alkali-metal atoms. Such improvement may 
possibly be achieved by augmenting powerful all-order 
techniques by contributions missed in a given order of 
many-body perturbation theory. We derived and ana- 
lyzed the entire set of fourth-order many-body diagrams 
for a one-particle operator. 

We highlighted the fourth-order contributions omitted 
in the popular coupled-cluster approach truncated at sin- 
gle and double excitations (CCSD). To recover the full 
set of fourth-order diagrams one should additionally con- 
sider the effect of triple excitations. In addition, the lin- 
earized version of CCSD should be augmented by non- 
linear contributions of double excitations. We presented 
explicit formulas for such complementary contributions 
in the Appendix. The representative diagrams may be 
found in Fig. ||. We also proposed a possible extension 
of the derived fourth-order contributions to all orders of 
perturbation theory. 

The derived expressions may be useful for an analysis 
of the completeness of all-order methods in the fourth 
order of perturbation theory and for designs of next- 
generation approximations in atomic many-body calcu- 
lations. 
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APPENDIX A: FOURTH-ORDER 
CORRECTIONS TO MATRIX ELEMENTS 

Here we tabulate fourth-order corrections to matrix 
elements of one-particle operator involving triple exci- 
tations and nonlinear contribution from double excita- 
tions. The classification of the diagrams and notation 
were introduced in the main text of the paper. Briefly, 
the matrix elements gijik of the Coulomb interaction are 
defined by Eq. (^. The quantities ijijik are antisymmet- 
ric combinations gijik — Qijik — gijki- Matrix elements of 
a non-scalar one-particle operator Z are denoted as z^ . 
Core orbitals are enumerated by letters a,b,c, d, com- 
plementary excited states are labelled by m,n,r,s, and 
valence orbitals are denoted by v and w. The notation 



are to be calculated by taking the hermitian conjugate of 
all preceding contributions and swapping labels v and w. 

For convenience of drawing the graphs, the sequence 
of interactions in numerators is sorted so that the inter- 
action to the right of another interaction appears lower 
in the corresponding Brueckner-Goldstone diagram. 



Zlx2{Tc) — 

X ^ 9abnr^cv9nrcm9mwab 



abcmnr 

E 

abcmnr 



+ 



Qabnr^cvQrwcmQmnab 
Qabiiv ^crdrwcmSrnnab 

abcmnr 



E 

abcmnr 

E 

abcmnr 

E 

abcmnr 



E 

abcmnr 



Qabrv ^cnQnrcmQmwab 



i^^mw 


^ab) {^rv ^ab) {^nrw 
Qabrv '^cnQrwcmQmnab 


^abc) 


(^mn 


^ab) i^^rv ^ab) i^^nrw 
Qbcnr^avQnrcmQmwab 


£abc) 


i^^mw 


Qbcnr^avQrwcmQmnab 


^abc) 


i^^mn 


^ab') i^^nr ^6c) i^^nrw 
Qbcnv^arQrwcmQmnab 


^a6c) 



, \^mn ^ab) \^nv ^bc) \^nrw ^abc) 

abcmnr 



E 

abcmnr 

E 

abcmnr 

E 

abcdmn 



bcrv ^anQnrcmQmwab 



{^mw ^ab) i^^rv S^bc) {^nrw £abc) 

Qbcrv^anQrwcmQmnab ^ 
{^mn ~ ^ab) i^^rv ^bc) i^^nrw ~ ^abc) 

dbdmn^cvQancdQmwab 

mw - Sab) (e mnw ~ £bcd) 

Qbdmn^cvQawcdQmnab 



- 

, , \^mn ^ab) \^mn ^bd) [^mnw ^bcd) 
abcdmn ' 



9bdmv ■^cnQancdQmwab 



E 

abcdmn 



Qbdnv ^cmQancdQmwab 



{^7nw ^ab) (^ni? ^bd) i^^mnw ^bcd) 
Qbdnv ^cmQawcdQmnab 



(^mn - Sab) {Snv - Sbd) [s. 

abcdmn 



+ 



, , \Smv Sbd) \Smw Sab) [Smnw Sbcd) 
abcdmn 



mnw Sbcd) 

Q cdmn Zbv Q ancdQmwab 

, , [SfYin Scd) [Suiw Sab) yS-mnw Sl^cd) 
abcdmn ' 



stands for 



■ Ez- The terms denoted c.c. 



E 

abcdmn 

E 

abcdmn 



Qcdmn^bvQawcddmnab 



{Smn Sab) (Smn Scd) {Smnw Sjjcd) 

9cdmv^bn9ancd9mwab 
{Smv ~ Scd) {Sjfiw Sab) (Smnw Sjycd) 



9 



- E 

abcdmn 



Qcdnv^bmQancdffmwab 



abcdmn 

- E 

abcmnr 



Qcdnv^hmQawcdQmnab 



^bcd) 

+ C.C. 



9bcmn ^wr9arbc9mnav 
(^mn ^av) (^mn ^6c) (^mnr ^bcv^ 

Qbcmn^wrdarcv ffmnab 

abcmnr 



+ 



+ 



E 

abcmnr 

E 

abcmnr 



Qbcmr ^rnQanbcQmwav 



Qbcmr^rnQancvQmwab 
{^mr ^6c) (s^miy ^ab) (s^mnn; ^bcv) 



+ 



+ 



- E 

abcmnr 

E 

abcmnr 

- 



Qbcnr^rmQa7ihcf]nvwav 



{^^mw ^av) (^nr ^6c) {^^mnw ^bcv) 



Qbcnr^rmQawbcQmnav 



ipmn ^av) i^^nr ^bc) i^^mnw ^bcv^ 

9bcnr ^rmdancv 9mwab 

, [^mw ^ab) \^nr ^bc) y^mnw ^6c 

abcmnr 



+ 



+ 



E 

abcmnr 

- E 



9bcnr^rm.9awcv9mnab 



{,^mn ^ab) {^nr ^6c) (^mnu; ^bcv^ 

Qbcnr ^wmf} arbc9 mnav 



+ 



abcmnr 



9bcnr^wm9arcv9mnab 



, l^mn ^a6j l^nr ^6cj (^mnr ^bcv) 

abcmnr 



+ 



+ 



- E 

abcdmn 

E 

abcdmn 



9bdmn^cd9cincv9Tnwab 



9bdmn^cd9awcv9mnab 



+ 



abcmnr 

E 

abcmnr 

E 
E 

abcdmn 

- E 

abcdmn 



abcdmn 



{^mn 


9bwmn^cr9o.rcv9fTinab 


^bcv) 


i^^mn 


^ab) (^mn ^bw) i^^mnr 
9bwnr^cm9arcv9mnab 


^bcv) 


i^^mn 


^ab) {^nr ^bw) {^mnr ~ 

9 cdrn n ^bd 9anbc9rnwav 


~ ^bcv) 




^cd) i^^mw ^av) i^^mnw 
9cdmn^bd9ancv9m'wab 


^bcv) 



{^mn ^cd) {^mw - Sab) {S mnw ^bcv) 



+ 



+ 



+ 



+ 



9cdmn^bd9o'Wbc9mnav 



{^mn ^ov) {^mn ^cd) {^mnw ^bcv^ 
9cdm.n ^bd9aw cv Qmnab 
abcdmn ^^"^^ ^afe) (^mn ^cd) {^mnw ^bcv^ 
9cwmn^br9arbc9mnav 



+ 



+ 



. j2 

, \^mn ^av) y^mn ^cw) \^mnr ^bcv) 

abcmnr 



- E 



9cwmn^br9arcv9mnab 



abcmnr 



■ E 

abcmnr 

■ E 

abcmnr 



^1x2(7?) 



(^mn ^ab) i^^mn ^cw^ {^mnr ^bcv^ 



+ 



9cwnr^bm9arbc9n 



{^mn ^av) (^nr ^cw) {^mnr ^bcv^ 

9c w n r b m 9a r c v9mnab 
{^mn ^ab) (^nr ^cw) {^mnr ^bcv^ 



+ 



9abns^sr9rwbm9mnav 



, \^mn ^av) \^ns ^ab) \^nrw ^abv) 

abmnrs 



E 

abmnrs 



9abnsZsr9rwmv9mnab 



{^mn ^afe) (^ns ^ab) i^^nrw ^abv^ 



+ 



+ 



E 

abmnrs 

E 

abmnrs 



9 a bns r9rs bm 9n 



- E 

abmni 

- E 

abmnrs 

E 

'.bmni 

- E 



abmnrs 



{^mn 


^a?;) (^ns ^ab} {^nrs 
9abns^wr9Tsmv9Tnnab 


^abv) 


(^mn 


^ab) (^ns ^ab^ {^nrs 

9abrs^sn9n,rbm9mwav 


~ ^abv) 


{^mw 


^av) i^^rs ^ab) {^nrw 
9abrs'^sn9nrmv9mwab 


^abv) 




9abrs^sn9rwbm9mnav 


^abv) 




9abrs^sn9rwmv9mnab 


^abv) 



+ 



(^mn ^ab) (^rs '^ab) (^^nrto ^abv^ 



+ 



+ 



+ 



+ 



abmnrs 

9abrs ^wn9rsbm9mnav 

, (^mn ^av) (,^rs ^afe) (^nrs ^abv) 

abmnrs 



E 

abmnrs 



9abrs^wn9rsmv9'mnab 



+ 



+ 



■ E 

abcmnr 



9acnr ^bc9nrbm9mwav 



9acnr^bc9rwbm9Tnnav 

, l,^mn ^avJ V^nr ^ac) [^nrw ^abv) 

abcmnr 



+ 



+ 



E 

abmnrs 



9 awn s^br9rs bm 9mnav 



+ 



,9 o, u ! r s ^bti 9t s hm, 9 m n a v 

, \^mn ^av) \^rs ^aw) \^nrs ^abv) 

abmnrs 



+ 



E 

abcmnr 

E 

abcmnr 

E 

abcmnr 

E 

abcmnr 



9bcnr^ac9nrbm9mwav 


{^mw 


^av) {^nr ^bc) {^nrw 
9bcnr ^ac9nrmv9'mwab 


^abv) 


{^mw 


^ab) (^nr ^6c) {,^nrw 
9bcnr ^ac9rwbm9mnav 


^abv) 




9bcnr^ac9rwmv9mnab 


^abv) 



(^mn ^ab) (^nr ^6c) (^nriu ■^abv) 



+ 



+ 



+ 



10 



- E 



Qhwns^arQrshmS'mnav 



abmnrs 



E 

abmnrs 



Qbwns ^ardrsmvdmnab 



+ 



E 

abmnrs 

E 

abmnrs 



Qhwr s^anQr sbmQmnav 
{^mn ^av) (^rs ^hw) {^nrs ^abv^ 

Qhwrs ^anQrsrn v !J rrinah 
{^mn ^afe) (^rs ^bw) i^^nrs ^abv) 



ex. 



Zoxs{Sv[Tv]) = 



E 



^wnQabrsQrsbmQmnav 





^v) {^mn ^av) {^nrs 


^abv^ 




"^wnQabrsQrsrnvQmnab 






£u) (^mn ^ab) {^nrs 


^abv} 




^wnQbcmr Q arhcQmnav 






^v) i^^mn ^av) {^mnr ~ 


~ ^bcv^ 




^wn9bcmr9arcv9mnab 





abmnrs 

E 

abmnrs 

E 

abcmnr 

, \^n ^v) \^mn ^ab) y^mnr ^bcv) 

abcmnr 



+ 



+ 



+ 



+ 



- E 

abcmnr 

- E 

abcmnr 



^wrQbcmnQarbcQmnav 



{Sf ^v) {^mn ^ov) (^mnr ^bcv^ 

^wr^bcmnQarcvQmnab 

^ws9abnr9rsbm9mnav 



+ 



+ 



abmnrs (^mn ^av) (^nrs ^a6t?) 

abmnrs ^^"^ (^mn ^ab) (^nrs ^abv) 



+ 



+ c.c. 



Zox3(£'„[7;]) = 
- E 



^bm9cdnv9ancd9mwab 



+ 



abcdmn 

- ^ 

, \^mn ^ab) i^nw ^bv) \^nrw ^abc) 

abcmnr 



^bn9acrv9rwcm9mnab 



+ 



E 

abcdmn 

- 



^bn9cdmv9cmcd9m.wab 



{^mw ^ab) {^nw i^^mnw ^bcd) 

^bn9cdmv9awcd9mnab 



+ 



, , \^mn ^ab) \^nw ^bv ) \^mnw ^bcd) 
abcdrnn 



Zbr9a'Cnv9n,rcm9mwab 



abcmnr ^^"^''^ ^afe) (.^rw {,^nrw ^abc) 



+ 



'^hr!JacML-!Jr,vcin!J,tuuib 



abcmnr ^^^^ ^ab) {^rw ^6v) {^nrw ^abc) 



E 

abcmnr 



^ c n 9abrv9rw cm 9 m 7iab 



{^mn ^ab) 



nw ^cv I V^nrto 



+ 



^cr9abnv9nrcm9mwab 



abcmnr ^^'"''"^ ^ab) {^rw ^cv^ {^nrw ^a6c) 



■ E 

abcmnr 



^cr9abnv9rwcm9mnab 



{^mn ^ab) {^rw ^cv) {^nrw ^a6c) 



^dm9bcnv9ancd9mwab 



+ 



+ 



abcdmn ^^^^ ^afc) (s^miu ^dv^ {^mnw ^fecd) 



+ 



Zdn9bcmv9o,ncd9mwab 



, , v^mzu ^ab) \^nw ^dv) \^mnw ^bcd) 
abcdmn ^ ^ 



- E 



^dn9bcmv9awcd9mnab 



abcdmn 



i^mn Safe) (^nto fi^dw) {^mnw ^bcd) 



+ C.C. 



Zox3{Sc[Tc\) = 

\ ^ ^bvSacnrQnrcmQmwab 



r-^ isw - £b) (e mi(j £ab) {^nrw ^abc^ 



+ 



E 

abcmnr 



- E 
E 

abcdmn 



^bv 9cicnr9rwcm9Tnnab 



{Syj £5) (s'mn ^ab) {^nrw ^abc) 

^bv9cdmn9ancd9mwab 
{^w ^b) {^mw ^ab) {^mnw ^bcd) 



+ 



+ 



E 

abcmnr 

E 

abcmnr 



^bv 9cd'm n 9 a wed 9 m nab 




^b") (^mn ^ab) i^^mnw 


- ^bcd) 




^cv9abnr9nrcm9mwab 






£c) {^mw ^ab) i^^nrw 


^abc) 




^cv9ctbnr9rwcm9Tnnab 





{Syj £^c) (^mn ^ab) {^^nrw ^abc) 



+ 



+ 



^dv9bcmn9ancd9mwab 



abcdmn {^mw ^ab) {^mnw ^6cd) 

^ ^ '^dv9bcmn9awcd9mnab 

, , ^d) \^mn ^ab) \^mnw ^bcd) 

abcdmn ^ 



+ 



c.c. 



ZoUD,[T^]) = 



^bn9cwmr9o,rcv9mnab 



/-^ {Smn - £ab) {Snw - £bv) {e mnr ^bcv^ 



- E 

abcmnr 

E 

abcmnr 

- E 

abcmnr 

- E 

abcmnr 

- E 

abcmnr 



Zbr9crmn9cincv9'm'wab 



i^mw ^ab) i^^rw ^6u) {^mnw ^bcv^ 



+ 



'^br9crmn9awcv9mnab 



+ 



{^mn ^ab) {^rw ' mnw ^bc 



^br9 c w m, n9'J^^r c 1 ^ 9 m n a b 
{,^mn ^ab) i^^rw (^mnr ^bcv^ 

^cn9bwmr9arbc9mnav 
{^mn ^av^ {^nw ^cv) {^mnr ^6cv) 

^cn9b'wmr9o,rcv9'mnab 
{^mn ^ab) {^nw ^cv) {^mnr ^bcv^ 

^cr9brmn9anbc9mwav 



+ 



+ 



abcmnr ^^^^ ^af) {^rw ^cv^ {^mnw ^bcv^ 

^ ^ Zcr9br7nn9ancv9m'wab 

, \^mw ^ab) \^rw ^cv) \^mnw ^bcv) 
abcmnr 



+ 



11 



- E 



^crfjbrmnfjawbcffmnav 



abcmnr 



■ E 

abcmnr 



{^mn ^af) {^rw ^cv) i^^mnw ^bcv^ 

^cr9brmn9awcv9mnab 
{^mn ^ab) {^rw ^cv) (^mnio ^bcv^ 



+ 



+ 



^bs 9 awnr Qr smv fjmnab 



E 

abcmnr 

E 

abcmnr 



^cr9bwmn9arbc9mnav 



Zcr9bwmn9arcv9'mnab 
{,^mn ^ab) (^rw ~ ^cv) {^mnr ^bcv^ 

^dm9bcdn9anbc9rnwav 



+ 



mnw ^bcv) 



+ 



- E 

abcdmn 



^dm9bcdn 9cincv 9Tn wab 



{^mw ^ab^ {^mw ^dv^ {^mnw ^bc 

^dn9bcdni9onbc9mwav 



+ 



, , x^mw ^av) \^'iiw ^dv) \^mnw ^bcv ) 
abcdmn 



E 

abcdmn 



^dn9bcdm9o,wbc9'fnnav 



+ 



E 

abcdmn 

- E 



'-dnJJbc<lni!JanL:c!Jniu:ab 



^dn9bcdm9awcv9'mnab 



abcdmn 



{^mn ^ab) {^nw ^dv^ {^mnw ^6cu) 



+ ex. 



E 



^an9bwrs9rsbm9mnav 



^as9bsnr9nrbm9mwav 

, l^mw ^avj V^siu ^av) \^nrw ^abv) 

abmnrs 



abmnrs 



- E 



^as9bsnr9T'wbm9'mnav 



abmnrs 



E 

abmnrs 



^as9b'wnr9rsbm9mnav ^ 
Zbn9o:Wrs9Tsbm9fnnav 



+ 



+ 



, \^mn ^av) [^nw ^bv) \^nrs ^abv) 

abmnrs 



+ 



abmnrs 



abmnrs 



- E 

abn 

E 

abmnrs 

E 

bmm 

E 

o6mnrs 

E 

abmnrs 

- E 

abmnrs 



^bn9awrs9rsmv9mnab 



{^mn ^ab) {^^nw ^bv) {^nrs ^abv) 



+ 



^bs9o.snr9nrb7n9mwav 


{^mw 


^av^ i^^sw ^bv) i^^nrw 
^bs 9 asnr 9nrmv9m,wab 


^abv) 


i^^mw 


£ab^ (^£sw ^bv} {^nrw 
^bs9cbsnr9Twbm9Tnnav 


^abv) 


{^mn 


^av^ {^sw ^bv^ {^nrw 

^bs9o,snr9rw7nv9mnab 


^abv^ 


{^mn 


^bs9a'wnr9rsbm9mnav 


^abv^ 



+ 



+ 



+ 



+ 



+ 



abcmnr 



abcmnr 



{^mn 




^abv) 




^cn9abcr9rwbm9mnav 






^av^ {^nw ^cv) i^nrw 






^cn9abcr9rw7ri v9rtinab 





{^mn S^afc) {^nw ^cv) {^nrw ^abv^ 



- E 

abmnrs 

E 

,bcm7 

E 

he 

E 

hcmr 

E 

,bcm7 

E 

abcmnr 

E 

abcmnr 



9abmr^cn9nrcv9mwab 



+ 



+ 



abcmnr 



abcmnr 



^cr9abcn9nrbm9mwav 


{^mw 


^af) (,^rw ^cv) i^^nrw 
^cr9abcn9nrmv9mwab 


^abv) 


{^mw 


^cr9cibcn9rwbm9'mnav 


^abv) 


{^mn 


^av) {^rw ^cv) (^nrxu 

■^cr9abcn9rwmv9mnab 


^abv^ 


{^mn 




^abv^ 



+ 



+ 



c.c. 



abcmnr 



+ 



- E 



9 a bnr ^cm9nrcv9m wab 



{^mw ^ab) {^nr ^ab) (^nr ^cu) 



+ 



abcmnr 

X ^ 9abnr^cm9rwcv9rrtnab 



, l^mn ^ab) \^nr ^ab) \^rw ^cv) 

abcmnr 



E 



abc7nnr 



9acmn^br9rwbc9n 



+ 



E 



fJacnir''bn!Jnj-bcfJnuvac 

(^mr ^ac) {^mw ^av^ (^nr ^6c) 



a6cmnr 

9acnr^bm9'^Tbc9mwav 

t y^mw ^av) \^nr ^acj l^^nr ^6cJ 

abcmnr 



9acnr^bm9rwbc9mnav 



, [^mn £av) i^nr ^ac) \^rw ^bc) 

abcmnr 



+ C.C. 



^Oxsi^nl) 

E 



■^an9bcmr9rwbc9mnav 



abcmnr 



- E 



Zar9bcmn9nrbc9fn'wa^! 



+ 



abcmnr 

V ^ ^ar9bcmn9rwbc9mnav 



abcmnr ^av) (^rtu ^^aij) (^^riu ■^fec) 



+ 



+ 



E 

abcmnr 



^bm 9o,cnr9nrcv 9Tnwa b 



{^mw ^ab) {^mw ^6v) (^nr ^cv^ 



+ 



^bn9acmr9'r'wcv9mnab 



- z = 

, y^mn ^ab) \^nw ^bv) \^rw ^cv) 

abcmnr 



+ 



E 

abcmnr 



^br9acmn9nrcv9mwab 



i^^mw ^ab^ i^^rw ^bv^ {^nr ^cu) 



12 



E 

abcmnr 



ZbrQ. 



acmnyrwcvQmnai 



C.C. 



Finally, the normalization correction due to valence 
triple excitations is defined as 

The correction to normalization may be represented as 
(the terms denoted c.c. are to be calculated by taking the 
hcrmitian conjugate of all preceding contributions) 



Qahnr 9nrmv 9mvab 



^ ^ Qbcmn Qancv Qmvah 

abcmn (^^^ ~" ^bc) (^7ni; ~^ ^ab) 



Qabnr Qnrbm Qmvav 

abmnr (^"^ ~ ^'i) ^^^^ ~ ^^fc) 
^ ^ Qbcmn 9anbc Qmvav 

abcmn 



abmnr 



QabnrQrvbmQmnav 
i^^mn ~ ^au) (^^nr ~ ^ab) 
QabnrQrvmvQmnab 



E 

abcmn 

E 

abcmn 



gbc 



lavbcQn 



- SbcY 



9b. 



'cmnHavcvHmnab 



(^mn ^ab) (^rj 



+ 



+ c.c. 



[1] C. S. Wood, S. C. Bennett, D. Cho, B. P. Masterson, [9] 
J. L. Roberts, C. E. Tanner, and C. E. Wieman, Science [10] 
275, 1759 (1997). 

[2] S. C. Bennett and C. E. Wieman, Phys. Rev. Lett. 82, [11] 
2484 (1999). 

[3] A. Derevianko, Phys. Rev. Lett. 85, 1618 (2000); V. A. [12] 
Dzuba et al, Phys. Rev. A 63, 044103 (2001); M. G. 
Kozlov et al, Phys. Rev. Lett. 86, 3260 (2001); W. R. [13] 
Johnson et al, Phys. Rev. Lett. 87, 233001 (2001); A . I. 
Milstein and O. P. Sushkov, e-print: ^iep-ph/0109257| ; A. [14] 
Derevianko. Phvs. Rev A. 6 5. 012106 (2002); Dzuba et 
al, e-print: |hep-ph/0111019t [15] 

[4] V. A. Dzuba, V. V. Flambaum, and O. P. Sushkov, Phys. 

Lett. A 141, 147 (1989). [16] 

[5] S. A. Blundell, W. R. Johnson, and J. Sapirstein, Phys. 
Rev. Lett. 65, 1411 (1990), Phys. Rev. D 45, 1602 (1992). 

[6] M. Y. Amusia and N. A. Cherepkov, Case Studies in [17] 
Atomic Physics 5, 47 (1975). 

[7] J. Sapirstein, Rev. Mod. Phys. 70, 55 (1998). [18] 

[8] W. R. Johnson, Z. W. Liu, and J. Sapirstein, At. Data 
Nucl. Data Tables 64, 279 (1996). 



J. Cizek, J. Chem. Phys. 45, 4256 (1966). 

F. Coester and H. G. Kiimmel, Nucl. Phys. 17, 477 

(1960). 

I. Lindgren and J. Morrison, Atomic Many-Body Theory 

(Springer- Verlag, Berlin, 1986), 2nd ed. 

S. A. BlundeU, W. R. Johnson, Z. W. Liu, and 

J. Sapirstein, Phys. Rev. A 40, 2233 (1989). 

S. A. Blundell, W. R. Johnson, and J. Sapirstein, Phys. 

Rev. A 43, 3407 (1991). 

M. S. Safronova, W. R. Johnson, and A. Derevianko, 
Phys. Rev. A 60, 4476 (1999). 

S. A. Blundell, D. S. Guo, W. R. Johnson, and 
J. Sapirstein, At. Data Nucl. Data Tables 37, 103 (1987). 
S. Wolfram, The Mathematica Book (Wolfram Me- 
dia/Cambridge University Press, Champaign, Illinois, 
1999), 4th ed. 

R. F. Bishop and H. G. Kiimmel, Physics Today 3, 52 
(1987). 

M. S. Safronova, A. Derevianko, and W. R. Johnson, 
Phys. Rev. A 58, 1016 (1998). 



